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Abstract 

Recent  developments  in  magnetic  switching  technology  have  enabled 
linear  induction  accelerators  to  operate  at  high  repetition  rates.  These 
accelerators  provide  high-average-power  electron  beams  for  many 
applications.  However,  some  applications,  such  as  free-electron  lasers,  also 
require  a  beam  with  a  very  narrow  energy  spread.  As  a  result,  the 
characteristics  of  the  waveform  driving  the  induction  cell  are  very  critical. 
Several  accelerator  systems  at  Lawrence  Livermore  National  Laboratory  use 
three-stage  magnedc  pulse  compressors  in  which  the  final  stage  consists  of  a 
pulse-forming  line  (PFL)  followed  by  a  saturable  output  reactor.  Ideally,  the 
PPL  forms  a  pulse  with  the  required  rise  time  and  flatness.  Unfortunately, 
traveling  waves  introduced  on  the  PFL  during  the  charging  process  as  well 
as  output  switch  performance  may  degrade  the  pulse  shape. 

To  understand  the  effect  of  these  factors  upon  the  output  pulse  shape, 
we  have  used  the  computer  code  SPICE  to  model  the  final  stage.  The  model 
simulates  PFL  dynamics  under  various  conditions  and  includes  a  first-order 
model  of  the  saturable  output  reactor.  Discussed  in  this  paper  are  the 
saturable  inductor  model  and  results  obtained  with  the  final  stage 
compression  model. 


Introduction 

Linear  induction  accelerators  (LIAs)  produce  intense,  high  energy, 
pulsed  electron  beams  for  use  in  many  applications.1-3  Typical  durations  of 
these  electron  beam  pulses  are  60-80  ns.  Such  accelerators  can  deliver  bursts 
of  electron  beam  pulses  occurring  at  multi-kilohertz  repetition  rates. 

An  LIA  accelerates  an  electron  beam  passing  through  a  series  of  ferrite- 
loaded  induction  cells.  A  voltage  pulse  driving  each  cell  induces  a  changing 
flux  within  the  ferrite  material.  The  changing  flux,  in  turn,  induces  an 
accelerating  potential  along  the  axis  of  the  electron  beam  and  increases  its 
energy. 

Induction  cell  drivers  employed  as  magnetic  pulse  compression 
networks4-6  enable  LIAs  to  operate  at  high  repetition  rates.  Unfortunately, 
these  networks  are  often  difficult  to  model,  but  models  are  essential  for 
network  design  and  optimization. 

We  have  used  the  SPICE  circuit  code7  to  simulate  the  dynamics  of  the 
last  stage  of  a  three-stage  magnetic  compression  network.  We  have 
developed  a  model  to  identify  and  optimize  parameters  important  to  the 
formation  of  the  accelerator  drive  pulse.  In  this  paper  we  describe  the 
magnetic  pulse  compressor,  the  SPICE  model,  and  results  we  obtained 
when  modeling  saturable-reactor  output  loss. 

Magnetic  Modulator  Description 

Figure  1  is  a  simplified  circuit  diagram  of  the  MAG-l-D  magnetic  pulse 
compressor  used  at  Lawrence  Livermore  National  Laboratory.  Capacitor  Cj 
is  initially  charged  to  25  kV  in  about  3.8  |xs,  at  which  time  L]  saturates  and 
allows  Cj  to  discharge  into  the  primary  of  the  1:10  step-up  transformer.  Lj 
has  been  designed  to  have  a  saturated  inductance  such  that  Ci  is  charged  to 
250  kV  in  about  1  pis. 

When  C2  is  fully  charged,  L2  saturates  and  begins  to  charge  the  2-12, 
35-ns,  water-filled  pulse-forming  line  (PFL).  The  charged  PFL  saturates 
L3,  whereupon  the  PFL  is  discharged  into  the  load.  The  modeled  portion  of 
the  network  is  enclosed  within  the  dotted  lines  of  Fig.  1.  Specifically,  the 
model  simulates  the  dynamics  of  the  PFL  charging  and  then  discharging  into 
the  load. 

The  output  switch  consists  of  a  single-turn  saturable  reactor  whose 
magnetic  material  consists  of  0.6-mil-thick  2605CO  Metglas.  The  flux 
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Figure  1.  Simplified  circuit  diagram  of  the  MAG-l-D  magnetic 
pulse  compressor. 


swing  of  this  material  is  approximately  3.2  T.  The  reactor  core  has  an 
effective  cross-sectional  area  and  mean  magnetic  path  length  of  60  cm2  and 
95  cm,  respectively. 

Three  physical  characteristics  of  the  output  reactor  must  be  represented 
as  quantities  suitable  for  the  SPICE  model.  Specifically,  these 
characteristics  are  the  unsaturated  inductance,  the  saturated  inductance,  and 
the  volt-time  product  of  the  reactor.  The  unsaturated  inductance,  L,^^,  is 
calculated  by  using  the  core  dimensions  and  estimating  the  relative 
permeability  of  the  core  based  on  the  magnetization  rate.  Similarly,  the 
saturated  inductance,  Lsat,  is  calculated  from  the  core  dimensions  and 
allowing  the  saturated  permeability  of  the  core  to  be  about  1.5  times  the 
permeability  of  free  space,  k0.  Finally,  the  volt-time  product,  <vt>,  of  the 
core  is  calculated  by 


<vt>  =  NAAB  , 


(1) 


where  N  is  the  number  of  turns  on  the  reactor,  A  is  the  effective  cross- 
sectional  area  of  the  reactor,  and  A B  represents  the  flux  swing  of  the  reactor 
core.  The  calculated  values  of  Lunsat,  Lsal,  and  <vt>  are  18.94  pH,  23.3 
nH,  and  22.06  mVs,  respectively. 

Circuit  Model  Description 

Figure  2  shows  the  modeled  first  compression  stage,  PFL,  and  output 
reactor.  Since  the  simulation  begins  with  the  discharge  of  C2  into  the  PFL, 


Figure  2.  SPICE  circuit  model  of  the  MAG-l-D  pulse  forming 
line  and  saturable  output  reactor. 


*  Work  performed  jointly  under  the  auspices  of  the  U.S.  Department  of  Energy  by  LLNL  under  W-7405-ENG-48  and  for  the  DOD  under  SDIO/ATC  MIPR 
No.  W31RPD-7-D4041. 
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we  represent  the  first-stage  compression  reactor  in  its  saturated  state.  In  the 
actual  circuit,  after  the  PFL  is  fully  charged,  L2  is  forward  saturated  to 
prevent  the  PFL  from  discharging  back  into  C2.  Therefore,  diode  Dj  has 
been  included  in  the  model  to  simulate  this  effect.  The  PFL  is  modeled  as  a 
2.1  -Q,  35-ns,  lossless  transmission  line. 

The  most  challenging  component  to  model  was  the  saturable  output 
reactor.  Unfortunately,  readily  available  versions  of  SPICE  did  not  contain  a 
nonlinear  magnetic  component  model.  However,  several  models  have  been 
developed  for  simulating  nonlinear  magnetic  components  for  switching 
power  supplies.8  These  models  form  the  basis  of  the  model  we  use  to 
simulate  output  reactor  operation. 

The  output  reactor  is  modeled  as  an  inductor,  L3,  in  parallel  with  a 
dependent  current  source,  Fj.  The  inductor  represents  the  unsaturated 
inductance  of  the  reactor.  Controlling  the  current  flow  through  F[  is  a  circuit 
that  simulates  the  saturation  dynamics  of  the  reactor.  The  control  circuit 
consists  of  an  RC  network  (R10,  C10)  that  integrates  the  voltage  across  the 
reactor  model.  The  integrated  voltage  is  proportional  to  the  reactor  flux  and 
is  subsequently  shaped  to  represent  the  saturation  dynamics  of  the  reactor. 
The  voltage  source-diode  combinations  Vsat(+)-D3  and  Vsat(_)-D4  simulate 
the  positive  and  negative  saturation  fluxes,  respectively.  As  long  as  the  time 
integral  of  the  voltage  is  between  Vsat(+j  and  Vsat(_j,  both  D3  and  D4  are 
reverse  biased,  preventing  any  current  flow  through  V7.  However,  when 
the  magnitude  of  the  integral  exceeds  either  Vsat(+j  or  Vsat(_),  current  flows 
through  V7  and  controls  current  flowing  through  Fj.  The  value  of  Rjj 
determines  the  current  flowing  through  V7,  and,  therefore,  determines  the 
saturated  inductance  of  the  reactor  model. 

We  define  several  constants  to  make  the  model  easier  to  implement. 
Before  the  voltage  across  the  reactor  model  is  integrated,  it  is  multiplied  by 
the  constant  Bj.  This  correction  is  necessary  because,  without  a 
multiplicative  constant,  the  magnitude  of  the  integral  would  be  on  the  order 
of  30  mVs,  which  would  produce  about  30  mV  at  node  11.  For  the  circuit 
to  function  as  required,  the  voltage  at  node  1 1  must  be  significantly  higher 
than  the  voltage  drop  across  diodes  D3  and  D4.  In  addition,  the  constant  fi2 
sets  the  ratio  between  the  current  flowing  through  V7  and  the  saturated 
current  flow  through  the  main  terminals  of  the  reactor  model. 

The  current  supplied  by  Fj  can  be  expressed  as 

4l(0=^fJvL3df-^Vsat.  (2) 

The  current  flowing  through  an  inductor  as  a  function  of  the  terminal 
voltage  is  given  by 

4i»=  l\VLdt  ■  (3) 


The  first  term  of  Eq.  (2)  can  be  rewritten  as 


dt  ; 


thus,  the  value  of  Rn  is  given  by 

Rn  =  BlB2Lsat  .  (4) 

In  this  model,  we  choose  Bj  and  B2  to  be  106  and  103,  respectively.  As 
a  result,  Ru  specifies  Lsat  in  nanohenries.  The  values  of  Vsat(+j  and  Vsat(.j 
are  related  to  the  saturation  flux  as 


^4at(+)  ^sat(-)  —  <vt> 

=  NAAB 

—  NA(Bsat(+^  —  2isat(„))  , 

where  Bsat(+)  and  Bsat(_j  are  the  positive  and  negative  values  of  Bsat  (in 
mVs),  respectively. 

An  initial-condition  voltage  on  the  integration  capacitor,  Ci0,  represents 
the  reactor  bias  point.  Typically,  this  inital  voltage  represents  half  the  reactor 
volt-time  product.  As  a  result,  the  reactor  is  seen  to  be  fully  reset  before  the 
simulation  begins. 


Dynamic  Loss  Model 

We  modeled  the  dynamic  loss  by  placing  a  voltage-dependent  current 
source,  Gj ,  across  the  terminals  of  the  reactor  model  and  by  making 


the  current  flowing  through  Gj  a  third-order  polynomial  function  of  the 
voltage  across  the  reactor  model  terminals.  Therefore,  a  nonlinear  resistor  is 
simulated.  When  a  voltage  is  present  across  the  reactor,  power  is  dissipated 
by  virtue  of  the  current  flowing  through  Gj.  Since  the  voltage  across  the 
reactor  is  directly  proportional  to  the  magnetization  rate,  dB/dt,  the  amount 
of  energy  dissipated  in  the  reactor  while  sustaining  a  voltage  is  also 
dependent  on  dB/dt. 

By  using  a  third-order  polynomial  to  specify  the  leakage  current,  we  can 
use  available  core  loss  data  to  simulate  reactor  loss.  Several  publications 
contain  loss  data  for  Metglas  2605CO9'10  for  magnetization  rates  of  up  to 
1 0  T/ps,  which  is,  unfortunately,  much  lower  than  the  magnetization  rate  of 
the  output  reactor.  Therefore,  we  estimate  losses  above  10  T/ps  to  be 
proportional  to  the  square  root  of  the  magnetization  rate.  The  loss,  once 
ascertained,  can  be  translated  into  a  leakage  current  as  a  function  of  the 
voltage  across  the  saturable  reactor.  We  then  use  a  curve-fitting  routine  to 
calculate  the  coefficients  of  the  polynomial.  For  the  loss  to  be  independent 
of  voltage  polarity,  the  polynomial  must  be  an  odd  function  of  voltage;  i.e., 
4:,i(v)  =  ~'gi(~v)- 

For  this  model,  we  set  the  coefficients  of  the  polynomial  at  0.0, 1.075  x 
10~2,  1.379  x  10-8  and  -1.46  x  10-7.  Figure  3  shows  the  leakage  current 
as  a  function  of  the  voltage  across  the  reactor  model  with  these  coefficients. 

Reactor  Model  Test  Circuit 

To  characterize  the  saturable  reactor  model,  we  modeled  the  test  circuit 
of  Fig.  4.  With  this  circuit,  the  saturated  inductance  and  the  loss  can  be 
calculated  and  compared  with  expected  results  before  the  model  is  used  with 
a  more  complex  circuit.  The  test  circuit  consists  of  a  capacitor  initially 
charged  to  a  voltage  V0.  Then,  capacitor  voltage  and  current  are  monitored 
(Fig.  5)  as  the  capacitor  is  discharged  through  the  reactor.  When  the 
capacitor  is  first  discharged,  current  flows  through  the  reactor,  setting  the 
core  to  +Bsat  and  reversing  the  voltage  on  the  capacitor.  The  reversed 
capacitor  voltage  appears  across  the  reactor  until  the  <vt>  of  the  reactor  is 
exceeded,  at  which  point  saturation  occurs  and  C]  begins  to  discharge 
through  the  reactor  once  again.  This  process  continues  until  the  energy 
initially  stored  in  Cx  has  been  dissipated. 

One  can  calculate  the  saturated  inductance  of  the  modeled  reactor  by 
measuring  the  base  width  of  the  current  pulse.  For  the  case  shown  in  Fig.  5, 
the  base  width  was  approximately  70  ns.  Since  the  test  circuit  is  a  simple, 
very  undamped  LCR  circuit,  the  saturated  inductance  is  easily  calculated  as 
24.3  nH.  In  addition,  values  for  energy  stored  in  the  capacitor  on  each 
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Figure  3.  Estimated  leakage  current  through  the  MAG-l-D 
output  reactor  as  a  function  of  voltage. 


Figure  4.  Test  circuit  used  to  measure  the  reactor  model. 
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Figure  5.  Voltage  and  current  waveforms  of  the  saturable 
reactor  model  in  the  test  circuit. 


successive  voltage  reversal  can  be  used  to  calculate  the  energy  loss  over  a 
range  of  dB/dt  values.  The  energy  loss  as  a  function  of  dB/dt  is  shown  in 
Fig.  6. 

The  test  circuit  can  also  be  used  to  produce  a  B-H  curve  by  plotting  the 
voltage  at  node  11  [V(ll)]  as  a  function  of  I(VCi).  V(ll)  is  directly 
proportional  to  the  reactor  flux,  B,  and  I(VCj)  is  directly  proportional  to  the 
magnetization  field,  H.  Scaling  V(ll)  by  the  inverse  of  the  effective  cross- 
sectional  area  of  the  core  yields  a  value  for  reactor  flux  in  Teslas.  Scaling 
ICVCj)  by  the  inverse  of  the  mean  magnetic  path  length  yields  a 
magnetization  field  value  in  units  of  ampere-turns  per  meter.  Figure  7 
shows  the  B-H  loop  we  calculated  for  the  test  reactor  model. 

Examples  of  Model  Results 

We  used  the  SPICE  model  of  the  last  section  of  the  pulse  compressor 
shown  in  Fig.  2  to  predict  the  characterisdcs  of  the  output  pulse  as  the  circuit 
parameters  are  varied.  Model  parameters  that  represent  the  present  pulse 
compressor  network  are  tabulated  in  Table  1. 


dB/dt  obtained  in  the  test  circuit. 


Capacitor  C2  is  initially  charged  to  250  kV.  The  results  of  the 
simulation  are  shown  in  Figs.  8  and  9.  Reflections  on  the  PEL  distort  the 
usual  (1  -  cosaw)  charging  waveform,  as  shown  in  Fig.  8.  These  reflections 
arise  because  the  PFL  charging  period  is  not  significantly  longer  than  the 
PFL  transit  time.  When  the  output  reactor  saturates,  the  PFL  is  discharged 
into  the  load,  as  shown  in  Fig.  9.  The  saturated  inductance  of  the  reactor 
model  limits  the  output  pulse  rise  time.  Charging  reflections  cause  the  ripple 
riding  on  top  of  the  pulse.  These  ripples  can  be  minimized  by  varying  L2 
slightly  (which  changes  the  PFL  charge  time)  and  discharging  the  PFL  at  the 
optimum  time. 

The  simulation  was  repeated  for  three  values  of  Lj  while  <vt>  was  held 
constant  at  22.06  mVs.  The  top  portions  of  the  output  pulses  are  compared 
in  Fig.  10.  As  can  be  seen,  output  flatness  is  a  strong  function  of  L2. 
Increasing  L2  from  625  nH  to  675  nH  increases  the  charge  time  from  236  ns 
to  247  ns.  Figure  1 1  shows  the  results  of  simulations  for  three  values  of 
<vt>  while  L2  was  held  constant  at  650  nH.  Once  again,  the  flatness  of  the 
output  pulse  is  strongly  influenced  by  the  precise  time  at  which  the  PFL 
discharges  into  the  load. 


Table  1.  Correlation  between  circuit  and  model  parameters. 


Circuit  parameters 

Model  parameters 

Component 

Value 

Component 

Value 

c2 

19.3  nF 

c2 

19.3  nF 

L2  (sat) 

650  nH 

l2 

650  nH 

Ti 

2.1  C2/35  ns 

Ti 

2.1  £2/35  ns 

L3  (unsat) 

18.94  pH 

l3 

18.94  pH 

L3  (sat) 

23.3  nil 

RU 

23.3  12 

L3  (Bsat) 

-1.6  T 

Vsat  (+) 

-11.03  V 

l3  (B0) 

1.6  T 

Veto  (t  =  0) 

11.03  V 

PFL  voltage 


Reactor  current  (kA) 

Figure  7.  A  B-H  curve  of  the  saturable  reactor  model.  The 
volt-time  axis  is  proportional  to  B.  The  reactor  current  axis  is 
proportional  to  H. 
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Figure  9.  PFL  output  pulse  across  2-12  load. 
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Conclusions 


Various  values  of  L2 
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Figure  10.  Magnitude  of  output  voltage  for  various  values  of 
the  PFL  charging  inductance  h2. 
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Figure  11.  Magnitude  of  output  voltage  for  various  values  of 
<vt>.  The  value  <vt>  is  proportional  to  AB  of  the  reactor  core. 


Various  values  of  (V+  ) 


We  have  used  the  SPICE  circuit  simulation  code  to  model  the 
characteristics  of  a  magnetically  switched  PEL.  We  modeled  the  output 
switch  with  a  subcircuit  that  simulated  the  output  reactor  saturation 
properties.  In  addition,  a  simple  loss  model  based  on  estimates  of  core  loss 
at  very  high  magnetization  rates  was  used  to  simulate  the  loss  characteristics 
of  the  reactor.  The  overall  model  has  been  useful  in  understanding  how  the 
shape  of  the  output  pulse  is  influenced  by  component  values  and  other 
parameters. 
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